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Abstract: We present results for the production of a top pair in association with a W-boson at next- 
to-leading order. We have implemented this process into the parton-levcl integrator MCFM including 
the decays of both the top quarks and the W-bosons with full spin correlations. Although the cross 
section for this process is small, it is a Standard Model source of same-sign lepton events that must 
be accounted for in many new physics searches. For a particular analysis of same-sign lepton events 
in which 6-quarks are also present, we investigate the effect of the NLO corrections as a function of 
the signal region cuts. 
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1 Introduction 

In this note we present results for the production of a top quark pair in association with a W- 
boson, calculated at next-to-lcading order (NLO) in QCD. The core of the calculation, performed at 
amplitude level, is essentially identical to the calculation performed for the bbW process in ref. [1], 
where results were given retaining the mass of the b quark. Since the calculation is performed at 
amplitude level we can treat the top quarks as on-shcll but include the effects of decay, retaining all 
spin correlations [2, 3]. The decay is included using a generalization of the amplitude method described 
in ref. [3]. The complete process has been implemented into MCFM v6.3. 

2 ttW production without decay 

To set the scene we display lowest order cross sections for several processes in which a tt pair is 
produced in association with a massive boson in Fig. 1. These are the total cross sections without 
branching ratios to observable modes and we consider LHC operating energies from 7 to 14 TeV. We 
note that, since we are presenting results for a proton-proton collider, the rates for tiW + and ttW~ 
are not equal. As well as the tt W process we have included the LO rates for tt Z and tt H, which are of 
a similar size. At the current energy of the LHC, y/s = 7,8 TeV, the cross sections are quite small, so 
that measurements of these processes may have to wait for data produced at the maximum energy of 
the LHC machine. Despite the fact that measurements have not yet been performed, these processes 
already play a role in searches for new physics. For this reason it is essential that the best possible 
theoretical predictions for these processes are made available. Next-to-leading order corrections to 
the Higgs signal process at the LHC were presented in Refs. [4-6] and the ttZ process has also been 
computed at this level [7-10]. Here we focus on results for the remaining ttW process. 

In passing we note that the production of other particles in association with a top quark pair 
has also been considered in the literature. For example, NLO results for the production of tt+jet 
were presented in ref. [11] and predictions for ttj were given in ref. [12]. The tibb process was also 
computed at NLO in refs. [13-16]. We shall not comment further on these processes since they are 
not the focus of this paper. 

Before proceeding to consider results for the tt W cross section including the decay, we present 
NLO results for the tt W cross section for a stable top quark and top antiquark. After adjusting the 
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Top pair associated production in lowest order 
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Figure 1. Lowest order cross sections for several processes in which a tt pair is produced in association with 
a boson at a pp collider, as a function of the centre of mass energy. Renormalization and factorization scales 
have been set equal to m t = 172.5 GeV and the ttH process is computed for a Higgs mass of 125 GeV. 

mass parameter, this calculation is essentially identical to our bbW calculation presented in ref. [1]. 
We note that a value for the NLO rate for the process tt W + at s/s = 7 TeV has been given in Ref. [17], 
although to the best of our knowledge no detailed phcnomenological studies have been performed. 
Throughout this paper we use the following electroweak parameters, 

M w = 80.398 GeV , IV = 2.1054 GeV , G F = 1.16639 x 10~ 5 GcV~ 2 , (2.1) 

and take the top and bottom quark pole masses to be, 

m t = 172.5 GeV , m b = 4.7 GeV . (2.2) 

For the parton distribution functions (pdfs) we use the sets of Martin, Stirling, Thorne and Watt [18]. 
For the calculation of the LO results presented here we employ the LO pdf fit, with 1-loop running 
of the strong coupling and ot s (Mz) = 0.13939. Similarly, at NLO we use the NLO pdf fit, with 
a s (M z ) = 0.12018 and 2-loop running. 

We begin by assessing the scale dependence of the predictions at LO and NLO. This is particularly 
interesting since one might imagine a rather wide possible choice of scales, say from mt to 2mt + mw- 
Our results for the LHC operating at s/s = 7, 8 and 14 TeV, (setting the renormalization scale /z r 
and factorization scale /i/ equal to a common scale /i), are shown in Figs. 2, 3 and 4 respectively. 
As expected, the scale dependence is reduced after including the NLO effects and the behavior of the 
predictions for ttW + and ttW~ is very similar. Qualitatively the scale dependence is also very similar 
at y/s = 7 and 8 TeV. However, at 14 TeV the contribution to the cross section from quark-gluon 
initial states, that enter the calculation for the first time at 0(a^), is much more important. Since 
this is a LO contribution the inherent scale dependence is uncancelled. As a result the improvement 
in scale dependence from LO to NLO is less dramatic at 14 TeV. 

We now present our best predictions for the cross sections for this process at LO and NLO, together 
with an assessment of the theoretical uncertainty at each order. Our predictions at the various LHC 
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Figure 2. Dependence of the total ttW cross sections on the renormalization and factorization scales jj, at 
— 7 TeV. The scale dependence is reduced going from LO (dashed curves) to NLO (solid curves). 




operating energies are given in Table 1. The central value of the cross sections corresponds to a scale 
fj, r = fif = mt and the central pdf set. The first uncertainty corresponds to a variation of the scales 
in the range [m t /4, 4toJ, which spans most of the conceivably relevant kinematic region. The second 
uncertainty originates from the pdf fitting procedure and is computed with the 90% confidence-level 
pdf sets and procedure of Ref. [18]. At NLO, where uncertainty sets are available that include both 
pdf and a s uncertainties, we use the method of Ref. [19] to also include the 90% confidence level 
uncertainty in a s (Mz) (corresponding to 0.1163 < a s (Mz) < 0.1233). Even considering the generous 
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Figure 4. As in Fig. 2 but at s/s = 14 TeV. 
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Table 1. Leading and NLO cross sections in femtobarns, with relative theoretical uncertainties estimated 
from scale dependence (in the range /j, — m t /4,m t ,4m ( ) and 90% confidence pdf uncertainty sets (including 
also ot s variation at NLO). 



90% confidence-level ranges we have considered, the uncertainty from the pdf and a s determination 
is smaller than the residual scale uncertainty of the NLO calculation. At 14 TeV this combined 
uncertainty is a factor of three smaller than the uncertainty from the scale dependence. In summary, 
although the cross section has a smaller uncertainty at NLO than at LO, the combined theoretical 
uncertainty is still sizeable (at best in the range of 30-40%) and grows with energy. 

3 ttW production including decay 

One of the advantages of the amplitude method is that it allows inclusion of the decay of the top quark 
and subsequent decay of the IF-boson, with all the correct spin correlations included. An example of 
a lowest order diagram, that displays the full decay chain that we include, is shown in Fig. 5. The 
inclusion of these decays is achieved with very little computational cost due to the fact that we can 
factorize all amplitudes into production and decay stages, see for example ref. [3]. In view of the 
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Figure 5. Example of a leading order diagram leading to a same-sign lepton event. The double lines indicate 
that the top and anti-top are on mass shell. 



smallness of the cross section we have not yet implemented radiative corrections in the decay of the 
top quark, although such an extension would be straightforward. 

Our calculation can also be used to assess the effects of radiative corrections in kinematic regions 
that are relevant for new physics searches. We shall study this in the following section, focussing in 
particular on same-sign lepton searches. 

4 Same-sign lepton events 

As can be seen from Fig. 5, the tiW process can lead to final states that contain same-sign leptons, a 
relatively rare phenomenon in the Standard Model. As such, these final states form the basis of many 
searches for theories of new physics, in particular supersymmetry, that could give rise to much higher 
rates. These searches can focus solely on the presence of two same-sign dileptons [20] or, more often, 
also require missing transverse momentum and jets [21-23]. The jets that are present may also be 
6-tagged [24]. 

Same-sign lepton events originating from the Standard Model are experimentally of three types. 
The first category is fakes, for instance when a hadron is misidentified as a lepton. The second, 
"q-flips" occurs when a lepton charge is mis-assigned so that an opposite-sign event enters the same- 
sign sample. Both of these categories are difficult to assess theoretically and must be estimated by 
data-driven techniques that rely on detailed detector simulation. The final category is comprised of 
rare Standard Model background processes that, in general, cannot be measured with current data 
samples. As a result, at present same-sign dilcpton analyses must rely on theoretical calculations for 
an assessment of this form of background. 

Examples of rare Standard Model backgrounds are: 

1. W + Z, W~ Z and ZZ production, with — > ^v, Z — > £~£ + . These processes can be computed 
at NLO accuracy including all spin correlations, see for instance Refs. [25, 26]. Predictions at 
NLO, interfaced with a parton shower framework, are also available [27, 28]. 

2. W^W^+dijets, with — > l^v. This process has been computed at NLO and included in the 
POWHEG framework [29]. 
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3. W^W*- production in vector boson fusion, with — > l^v. The NLO corrections to this 
process were presented in Rcf. [30]. 

4. WWW, W + W~Z, ZZZ in various decay channels, producing same-sign leptons with either 
additional leptons or jets. Predictions for these processes are available at NLO [31, 32]. 

5. ti Z production, with at least one of the top quarks decaying scmi-leptonically. This cross-section 
is available at the NLO level [7-10], but unlike our ti W calculation the complete spin correlations 
are not included at NLO. 

6. ti W production, with one top quark decaying scmi-lcptonically and the other hadronically. The 
prediction of ti W processes at NLO with subsequent decays (such as the one in the previous 
sentence) is the subject of this paper. 

The relative size of these various backgrounds is dictated by the further event selections that are 
applied in the search. If the signal region is defined by the presence of additional jets, for instance, 
then the diboson backgrounds shown above are relatively suppressed. Requiring that such jets are 
6-tagged favors both of the associated top production processes, where two 6-quarks are naturally 
produced. 




Figure 6. The distributions of fir (left) and Ht (right) for the tiW + process, with the top quark decaying 
leptonically and the antitop quark hadronically. The calculation is performed for the 8 TeV LHC, under the 
same-sign dilepton cuts described in the text. The LO prediction is the dashed (blue) histogram while the 
NLO prediction is solid (red). 



As an example that is especially sensitive to the ti W Standard Model process, we focus on a recent 
CMS analysis that requires the presence of same-sign dilcptons, missing transverse momentum and 
6-tagged jets [24]. We will provide predictions for the tiW background, with one top quark decaying 
leptonically and the other hadronically, that would be relevant for this search if it were repeated at 
8 TeV. Adopting the cuts from Rcf. [24], we require two same-sign leptons in the acceptance range, 

p T (lepton) > 20 GeV , |r/(lepton)| < 2.4 . (4.1) 
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Table 2. Definition of signal regions and cross sections in femtobarns for a single flavor of lepton in each 
lcptonic W decay. The first quoted uncertainty corresponds to variation of the scale by a factor of four in 
each direction about the central value, /i = m t . The second range corresponds to the 90% confidence level pdf 
uncertainty (LO) or pdf+a s uncertainty (NLO). 



Jets are clustered according to the anti-fey algorithm with a distance parameter D = 0.5 and we 
require that the algorithm finds at least two jets that satisfy, 

pr(jet) > 40 GeV , |r/(jet)| < 2.5 . (4.2) 

Two of these jets should contain the b and b quarks. Finally, we require that the missing transverse 
momentum, which in our calculation is given by the magnitude of the vector sum of the two neutrino 
four-momenta, is greater than 30 GeV. 

We follow the division of events into signal regions that is presented in Ref. [24] . These regions are 
defined by minimum cuts on the missing transverse momentum and the variable Ht that is defined 
by, 

H T = 5>T(jet i) , (4.3) 

i 

i.e. it is the scalar sum of the jet transverse momenta. Since these observables are crucial in the 
definition of the signal regions, in Fig. 6 we present our predictions for them obtained at LO and 
NLO. These predictions are for a pair of positively-charged leptons, i.e. the ttW + process, although 
they do not differ greatly for the negatively-charged pair (tiW~). We see that, although the shape of 
the missing transverse momentum distribution is unaffected by radiative corrections, the shape of the 
Ht distribution is changed considerably. At NLO the peak around 250 GeV is reduced but the tail is 
somewhat higher. 

Our results for the various signal regions at 8 TeV are presented in Table 2. Note that we show 
predictions for positively- and negatively-charged lepton pairs separately and not the sum. As such, 
regions SRI and SR2 of Ref. [24] are conflated. Moreover, our process contains only two 6-quarks and 
so does not contribution to the region SR7. The change in cross section from LO to NLO in each 
of the signal regions is relatively mild but the reduction in theoretical uncertainty from the choice of 
scale is substantial. Both the change in cross section and the associated scale uncertainty are largest 
for SR5, one of the regions defined by the largest Ht cut. Since the f,T cut is relatively small for this 
region, these results can be explained by the fact that the calculation is more susceptible to the change 
in shape of the Ht distribution shown in Fig. 6. We explore the effect of the NLO corrections to the 
positively-charged dilepton process more generally in Fig. 7. The figure shows the -ftf-factor, i.e. the 
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ratio of the NLO and LO predictions, as a function of the fir and Ht cuts in the range 30-130 GeV 
(fir) and 80-500 GeV (Ht)- The A"- factor is evaluated at our central scale choice, fi r = fif = m t . 
For most of these choices the A-factor is close to unity, with the value rising as the Ht cut is raised 
and the JSt cut is reduced. 
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Figure 7. The dependence of the A-factor for the positively-charged dilepton process on the cuts that are 
applied to fir and Ht- The A- factor is defined for /j, — mt, 



5 Conclusions 

We have presented first NLO results on the process tt including the decay of the top quark and 
the vector bosons. The cross section has interest, both as a Standard Model measurement, and as a 
source of events that contain same-sign lcptons, missing energy, jets and 6-quarks. 

At all potential LHC operating energies, between 7 and 8 TeV, evaluating cross sections with a 
central scale choice of [i = mt leads to very little enhancement of the LO cross section at NLO. For 
a large scale choice, such as ^ = 2m t + mw the A- factor is in the region 1.3-1.4. Considering scale 
variation and the combined pdf and a s uncertainty, the overall accuracy of the NLO prediction is at 
the level of 30% at best. 

At NLO the Ht observable, defined as the scalar sum of all jet transverse momenta in the event, is 
subject to important NLO corrections. These cause more events to be found at high Ht so that NLO 
effects arc more pronounced as progressively harder cuts on Ht are performed. This is of particular 
interest in the same-sign dilepton channel, where cuts of exactly this nature are usually performed in 
order to maximize the sensitivity to new physics. 
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